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Group theoretical methods are used to determine the electromagnetic properties of artificial mag- 
netic meta-materials, based solely upon the symmetries of the underlying constituent particles. 
Point groups for such materials are determined. From the transformation properties of an elec- 
tromagnetic (EM) basis under symmetries of the particles, it is possible to determine, (i) the EM 
modes of the particles, (ii) the form of constitutive relations (iii) magneto-optical response of a 
meta-material or lack thereof. These methods are shown to be useful for determination of the 
isotropic or bi-anisotropic nature of artificial magnetic particles. The results for several artificial 
magnetic metamaterials are given. Based upon these methods, we predict an ideal planar artificial 
meta-material, which eliminates an undesirable electric resonance while still exhibiting a magnetic 
response. Further we determine the subset of point groups of which particles must belong to in 
order to yield an isotropic 3D magnetic response, and we show an example. 



Progress in the rapidly advancing field of metamate- 
rials has experienced even further growth due to the re- 
cent observation of a negative index (NI) medium. j3| A 
NI medium is an artificial material in which both the 
magnetic fi and electric e response obtain simultaneous 
negative values, thus yielding a negative index of refrac- 
tion. Although predicted over three decades ago,0 it 
wasn't until recently that such materials were realized. 
Typically NI media utilize individual components for the 
electric and magnetic response. Materials which exhib- 
ited negative epsilon, known as artificial dielectrics, have 
been known since the 40's.|3, 0, The critical compo- 
nent responsible for the demonstration of NI was the re- 
alization that artificial non-magnetic materials could be 
constructed to exhibit negative response. 

The first demonstrations of NI materials were per- 
formed at microwave frequencies, due in part to the 
ease of fabrication as well as simplicity of measurements. 

The magnetic component of NI materials have since 
been demonstrated at lower RfQ frequencies and higher 
THzj^ frequencies. Great interest for the further exten- 
sion of these exotic materials to optical frequencies (uti- 
lizing nano-sized elements) is fuelling massive research 
efforts. Given the potential of NI materials to span the 
electromagnetic spectrum, it is important to understand 
their full complex electromagnetic behavior. In particu- 
lar, the structures utilized for NI are bianisotropic and 
may yield rich electromagnetic properties such as chiral- 
ity K and/or non-reciprocity x- 

At microwave frequencies where complex reflection and 
transmission measurements (S-parameter) are common, 
it is difficult to characterize bianisotropic materials. This 
difficulty stems from the fact that these materials are nec- 
essarily described by the most general form of the con- 
stitutive relations. There may be 36 complex quantities 
to determine, and thus standard complex reflection and 
transmission measurements yield incomplete EM infor- 
mation. At THz and higher frequencies where phase sen- 
sitive measurements are not common it is even more diffi- 



cult, although methods such as ellipsometry or THz time 
domain spectroscopy may prove useful. Further since a 
magnetic and electric resonant response enter similarly 
into the Fresnel equations, it is difficult to determine 
the origin of the response. Typically, artificial magnetic 
metamaterials are constructed from conducting elements, 
and thus a full electromagnetic characterization is a ne- 
cessity, since an electric response is unavoidable. An- 
alytical theories have yielded equations able to predict 
the resonant frequency and plasma frequency ujp of 
metamaterials. However there is a lack of suitable ana- 
lytical methods capable of determining the many varying 
and complicated EM properties. Simulation is still heav- 
ily relied upon and, as of yet, is still unable to determine 
and distinguishing bianisotropic response. 

We present a group theoretical method capable of de- 
termining various electromagnetic properties for mag- 
netic metamaterials. This method is based simply upon 
the symmetry operations of the constituent particles 
about a point in space, i.e. point group theory. An 
EM basis is assigned to the artificial particle under in- 
vestigation, and transformations of this basis under the 
symmetries of the group yield the electromagnetic modes. 
This new method allows one to calculate: the EM modes, 
the form of the constitutive relations, and the determine 
wether a particle will exhibit magneto-optical response. 
This approach is demonstrated by way of an example. 

Materials utilized for negative mag netic response have 
been shown to be bianisotropic, [lO| thus let us review 
the constitutive relations for these materials, which can 
be written as: [ill 
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In Eq^ the permittivities within the 2x2 matrix are 
tensors of second rank, also called dyadics. The terms 
^ and ^ are called the magneto-optical permittivities. 
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FIG. 1: Point group symmetries of the SRR particle. Panel 
(a) shows the coordinate system convention. Panel (b) shows 
the symmetry axis of the SRR and the C2 symmetry (rotation 
about the axis by 27r/n, n=2). Panels (c) and (d) demonstrate 
the mirror plane symmetries. 



and they describe coupling of the magnetic to electric 
response and electric to magnetic response respectively. 

As an example, let us start by examining the point 
group symmetries of a typical element utilized for mag- 
netic response. In Fig. ^ (b-d) the symmetry oper- 
ations of a split ring resonator (SRR) are shown. A 
symmetry operation brings the element into self coin- 
cidence. Thus for the SRR it can be seen that there 
are three symmetries which meet this criterion. In ad- 
dition, all groups also contain the identity operation E. 
We then find that the SRR particle belongs to the C2V 
point group which contains the following elements [E, C2, 
a{xy), <t{xz)] where we use the Schoenflies notation. 

We now turn to determination of the EM modes of the 
particles. Since we have determined SRRs to belong to 
the C2V group we can utilize the character table for this 
analysis(see Table I). With the character table we can 
assign a basis set to the SRR and see how this basis set 
transforms under the symmetry operations of the SRR. 
For brevity we only consider the outer split ring, how- 
ever the inner ring is also easily handled by this method. 
We want our basis to represent areas of electrical activ- 
ity, thus we choose the basis shown in the left column of 
Fig. Olb). Regions marked with arrows represent areas 
which can be polarized by an external electric field (i.e. 
currents can flow in these directions). These are similar 
to the P orbitals used in molecular orbital group theory 
(MOGT). The next step is to write out matrices which 



describe how this basis transforms under C2V For ex- 
ample, there are five vectors which make up our basis 
for the SRR, and since the identity leaves the particle 
unchanged, this would be a 5x5 matrix with I's along 
the diagonal. Once we obtain matrix representations for 
each element of the group, we can use a result of the Or- 
thogonality Theorem to determine how many times each 
irreducible representation (irrep) occurs. For our chosen 
basis we use the following equation, |0 

"™ ^ X! "-cX(5)Xm(5) (2) 

c 

where h is the order of the group (h=4 in this case), 
Uc is the number of symmetry operation in each class, 
xig) are the characters of the original representation, and 
Xm are the characters of the m*'' irreducible representa- 
tion. Using Equation (01 we find our basis is spanned by 
rsfl,fl=2Ai+3B2. 
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TABLE I: Character table for the C2v point group. The body 
of the table lists the characters, (trace of a matrix represen- 
tation) of the group. The first column lists the irreps of the 
group, the top of columns 2-5 lists the symmetry operations, 
and the last 2 columns lists a number of linear and quadratic 
functions that transform as the various irreps of the group. 
The bottom row lists the characters of our chosen SRR basis. 

Once we know the irreps spanned by an arbitrary basis 
set, we can work out the appropriate linear combinations 
of basis functions that transform the matrix representa- 
tions of our original representation into block diagonal 
form. These are called symmetry adapted linear combi- 
nations (SALCs). We use a projection operator to de- 
termine the SALCs, that transforms as an irrep. This is 
given by, 

<l>i ^^Xk{g)g(l)z (3) 

a 

where (f)'^ is the SALC, Xk is the character of the k*'' irrep, 
g is the symmetry operation, and (pi is the basis function. 

We can normalize the SALCs determined by Equa- 
tion ||2J), but this is not necessary since a constant factor 
doesn't affect the symmetry of the calculated modes. The 
response of the SRRs can now be determined by consider- 
ing incident external electromagnetic fields. For example 
by examination of the character table for the C2V point 
group we see that light polarized along the i-axis trans- 
forms as Ai, since it transforms in the same manner as 
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FIG. 2: Basis utilized (red arrows in left column) for magnetic metamaterials used to calculate the EM modes (red arrows in 
the remaining columns). The remaining columns for each row show the SALCs and modes of the SRR particle, as determined 
by point group symmetries. For each column, the irrep is shown above and the corresponding component of the an external 
electromagnetic wave, or function is shown. 



the function x. Thus y-polarized light transforms as B2 
symmetry. The function Rq, represents rotation about 
the a axis, where a=x, y, z. Thus a magnetic field polar- 
ized along the z-direction also transforms as B2 symme- 
try. We summarize these results in Fig. |21for: (a) planar 
spirals, (b) SRRs, (c) Omega particles 13] (d) an SRR 
and its enantiomer|l^ (e) symmetric ring resonator. 

An electric field polarized along the x-axis (Ai) of the 
SRR drives currents as shown in row (b) of Fig. |21 This 
would give a frequency response determined by the di- 
mensions (length) of the SRR segments along which Ej^ 
lies. For y-polarized (B2) light much more exotic behav- 
ior is predicted. Notice that for the irrep of B2 both y and 
R^ form a suitable basis. Thus we can use a linear combi- 
nation of these two functions for the basis. This predicts 
that the SRR should exhibit a magneto-optical response, 
as in accord with MOGT. 12] Furthermore, since Ej, and 
Bj: are the same basis, they will occur at the same fre- 



quency. In other words, an E field polarized along the 
y-axis will result in a resonant response at a frequency 
and a magnetic field polarized along the z-axis will 
result in a resonant response at the same frequency loq. 
These theoretical predictions are consistent with results 
obtained from a simple analytical model^^l as well as by 
simulation. |l5j 

In Fig. |2] we show predictions for other various mag- 
netic metamaterials. Particles with the lowest symmetry 
are listed on the top and higher symmetry particles on 
the bottom. The artificial magnetic metamaterials listed 
in the first 3 rows are shown to be bianisotropic, and thus 
we list the form of the constitutive equations governing 
the predicted magneto-optical response on the right side 
of each row. In row (d) we show a particular way of sym- 
metrizing the SRR particle, which we predict will elim- 
inate the magneto-optical response. The material is a 
bipartite lattice with each of the two sub-lattices consist- 
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ing of SRRs each with the gap oriented oppositely. Thus 
it is predicted that any polarization rotation or mixing 
resulting from one unit cell is corrected by the other unit 
cell, resulting in no net polarization rotation. Indeed it 
can be seen that the theory predicts no magneto-optical 
activity, while still exhibiting a magnetic response for B^. 
Another predicted way to eliminate ^ and C is to add a 
second gap in the SRR opposite to the first, as shown in 
Fig- Ete). Again the theory predicts no magneto-optical 
response and thus we have eliminated bianisotropy by 
symmetrizing the SRR. This structure is simpler than 
that depicted in Fig. Efd) and we have a predicted mag- 
netic response for B^. 

Lastly let us turn to the prediction of an ideal magnetic 
metamaterial. With our new understanding of the irreps 
of a point group and their relation to magneto-optical 
behavior we can leverage point group tables for help. An 
ideal magnetic particle is one in which no magneto opti- 
cal behavior is predicted. In the parlance of group the- 
ory this implies we should look for a point group which 
has linear basis functions with rotational functions R^ 
but with little or no occurrences of linear x, y, z func- 
tions. This not only ensures the elimination of the ^ and 
^ terms, but also off-diagonal terms in the /i and e re- 
sponse functions, like those which occur for the particle 
listed in Fig. E^d), i.e. 3Bg and 2B„. In Table II we 
show a candidate point group which should have good 
magnetic response with no magneto-optical activity, and 
no frequency dependent e occurring near the magnetic 
resonance. A particle which has the symmetry of this 
group is shown if Fig. O 

Group theoretical analysis is carried out for the par- 
ticle depicted in Fig. Ol^a), and we find the following 
modes r=Aig-|-A2g-l-Big-|-B2g-l-2E„. Thus the only lin- 
ear modes determined are a magnetic mode (A2g) and an 
electric mode (E„). The fact that the electric mode does 
not occur in the same irrep as the magnetic mode ensures 
it will not occur at the same frequency. Further the two 
dimensional E„ mode implies the electric response of the 
particle will be identical along the x and y directions with 
no cross coupling terms {exy=eyx=Q) ■ Thus if we want 
to construct a 3D isotropic magnetic metamaterial free 
from magneto-optical activity, then the geometry of the 
constituent particles are required to be one of the follow- 
ing point groups: Th,Td,^h, and Oh- The two simplest of 
which to visualize is I/j, an icosahedron, and O^, a cube 
(depicted in Fig. |3fb)) or octahedron. 

We have demonstrated a new method capable of deter- 
mining various properties of the electromagnetic response 
for artificial magnetic metamaterials. A specific example 
has been worked out for the most common element uti- 
lized for negative magnetic response. This analysis has 
also been carried out and detailed for other various arti- 
ficial magnetic metamaterials. We have predicted 2 ideal 
magnetic particles, including one for 3D isotropic mag- 
netic response. Whether one wants to take advantage of 
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TABLE II: Character table for the D4h point group. 




FIG. 3: Predicted ideal planar and ideal 3D magnetic parti- 
cles. In panel (a) we show a planar magnetic particle with 
D4h symmetry. The currents under the A2g magnetic mode 
are shown. Panel (b) shows a 3D isotropic magnetic particle 
with Oh symmetry. The electric response of this particles is 
also isotropic, but importantly does not occur at the same 
frequency as the magnetic resonance. 



the exotic electromagnetic properties that emerge from 
the bianisotropic nature of artificial metamaterials, or 
construct artificial materials free from this complication, 
these novel methods are valuable for determining the ex- 
pected response. 
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